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We present the conditions under which probabilistic error-free discrimination of mixed states is 
possible, and provide upper and lower bounds on the maximum probability of success for the case of 
two mixed states. We solve certain special cases exactly, and demonstrate how the problems of state 
filtering and state comparison can be recast as problems of mixed state unambiguous discrimination. 
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I. INTRODUCTION 



A characteristic feature of quantum mechanics is that 
if a system is prepared in one of a set of non-orthogonal 
pure states, then there is no measurement that can yield 
an error-free determination of which state was prepared. 
Nonetheless, it is possible to achieve a probabilistic error- 
free discrimination, that is, one which sometimes fails but 
when successful never gives an erroneous result. This 
sort of discrimination procedure is generally referred to 
as unambiguous discrimination (UD). The UD of two 
pure states prepared with equal prior probabilities was 
considered by Ivanovic and Dieks and the optimal 
procedure was given by Peres This was generalized 
to the case of unequal prior probabilities by Jaeger and 
Shimony_p| . The problem of three pure states was anal- 
ysed in 0| , while multiple pure states were considered in 
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It is a common misconception that the unambiguous 
discrimination of mixed states is impossible 0, 0- In- 
deed, in it is explicitly stated that "one cannot unam- 
biguously discriminate mixed states." That such a claim 
cannot be correct for an arbitrary mixed-state ensem- 
ble is proven, however, by the following counterexample: 
any set of orthogonal mixed states can always be discrim- 
inated with zero probability of error. 

What is less obvious is that there exist sets of non- 
orthogonal mixed states for which UD is possible. The 
critical feature of such sets is that their elements do not 
have identical supports |Ty| . In fact, all that is required 
for there to be a non-zero probability of error-free dis- 
crimination is that one of the density operators have a 
non-zero overlap with the intersection of the kernels of 
the others. 

In this paper, we consider the problem of determining 
the optimal UD procedure for an arbitrary pair of mixed 
states. We derive strong upper and lower bounds on the 
probability of a conclusive result, and we provide an exact 
solution in the special case where both states have kernels 
of dimension 1. 

The ability to unambiguously discriminate a pair of 
mixed states has many applications. Indeed, several re- 
cently studied problems can be recast as special cases of 



UD of mixed states. We shall focus on two tasks of this 
sort: state comparison [f| (determining whether two sys- 
tems are described by the same or different pure states) 
and state filtering pTol Hl| (discriminating one pure state 
from a set of pure states) . We find that for the state com- 
parison problem, our upper and lower bounds coincide, 
and thus automatically yield the optimal solution. For 
the state filtering problem, we find that our lower bound 
is equal to the optimal solution found in 0, • Given 
that our lower bound has a simple geometric interpreta- 
tion, it serves to clarify the nature of the optimal state 
filtering procedure. In particular, it shows that the eigen- 
bases of the optimal POVM elements depend only on the 
subspace spanned by the states against which one is try- 
ing to filter, and not on the specific states themselves. In 
addition, our approach generalizes in a straightforward 
way to more complicated estimation tasks, such as state 
filtering and state comparison when the unknown states 
are themselves mixed. 

It should be noted that our results apply not only to 
mixtures that arise from ignorance about which of several 
pure state descriptions applies, but also to those arising 
as the reduced density operator of an entangled state. As 
such, our results can be applied to the task of achieving 
an UD of two entangled states of a composite system 
given access to only one of the subsystems. It is also 
worth noting that our lower bound is obtained by making 
implicit use of the CS decomposition, which constitutes a 
powerful tool in both modern linear algebra and classical 
signal analysis To our knowledge, this is the first 

application of the CS decomposition in quantum signal 
analysis. 



II. GENERAL FORMULATION 

We consider the task of discriminating unambiguously 
between two mixed states po and pi with prior probabil- 
ities po and p\. The measurement procedure can have 
up to three outcomes, associated with identifying the 
state as po, identifying the state as pi, and failing to 
identify the state conclusively. The most general three- 
outcome measurement is represented in quantum me- 
chanics by a three-element positive-operator valued mea- 
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sure (POVM), which we denote by {E , Ei,E?}. Because 
the identification must never be in error, we require that 



Tr(po^i) = Tr{ Pl E ) = 0. 
The probability P of successful UD is: 

P = p Tr(p E ) +p 1 Tr(pxE 1 ). 



(1) 



(2) 



We shall denote the support of pb by supp(/9b) and its 
kernel by ICb (&=0,1). It is clear that any intersection of 
/Co and 1C\ is not useful for the purposes of discriminating 
Po and pi since neither state has any overlap with this 
subspace. We shall therefore assume, in what follows, 
that the Hilbert space is equal to the span of the supports 
of po and p\. A necessary and sufficient condition for 
satisfying is that the POVM element E (Ei) have 
support only in the subspace K.x{K.q). It follows that for 
there to be a non-zero probability of success, at least one 
of /Co an d K.\ must be non-zero. This occurs if and only 
if supp(po) ^ supp(pi). 

We seek to maximize subject to Q and the con- 
straint that Eq, Ei, and E? be positive and sum to the 
identity. It suffices to vary over positive Eq and E\ which 
satisfy 



I — E — Ei > 0. 



(3) 



This optimization problem is an instance of a semi- 
definite programming problem for which there exist effi- 
cient numerical algorithms. (See |6| for applications of 
semi-definite programming to UD of pure states, and 
0, @ f° r applications to the problem of optimizing the 
discrimination probability between mixed states for a 
fixed error rate, which interpolates between maximum 
likelihood estimation and UD). 



III. SOLUTIONS IN CERTAIN SPECIAL CASES 
A. Mixed states with orthogonal kernels 

If /C _L /Ci, the optimal POVM which satisfies © 
is clearly E = K\, E\ = Kq, where Kb is the pro- 
jector onto ICb- This yields a probability of success 
pmax = po Tr(p ^i) +piTr(pxK ), Note that this so- 
lution also applies when one of the kernels is the null 
space. Note also that commuting mixed states necessar- 
ily have orthogonal kernels. It follows that this result 
specifies the maximum probability of UD for overlapping 
classical probability distributions. 



B. Mixed states with 1-dimensional kernels 

We now turn to the special case wherein the kernels are 
both 1-dimensional - that is, the states po and p\ have 
rank n — 1, and the span of their supports is an n dimen- 
sional space. Denoting by |fc&)(fcf,| the projector onto ICb, 



the POVM must be of the form {E = a|fc 1 )(fci|,£' 1 = 
(3\k )(k \,E? = J-a|fci)(fci|-/?|fco)(fco|}. Our task there- 
fore becomes to compute 



P max = max 



nax (apo(ki\p \h) + Ppi{k \ P i\k ) ) , (4) 

where the maximization is subject to the constraint 

J-alAiXfcil-flAoXfcol >0. (5) 

For convenience, we define 9 to be the angle in Hilbert 
space between the 1-dimensional kernels, 



COS 2 EE |(fc |fcl)f . 



(6) 



Taking the minimum eigenvalue of the left hand side of 
iJSJ, we can re-express this inequality as 



cos z 



< 1. 



(7) 



Solving for p max under this constraint we find that 
a opt = (1 - y/A 1 /A Q cos 0)/ sin 2 9, (3 opt = (1 - 



yjA I A x cos 9) I sin 2 9, and 



P max = \ iT^e 11 CObU < V (8) 



where A = p Q (k 1 \p Q \k 1 ), A x = Pi<fco|pi|^o)j an d where 
A mia = min{A , Ai}, A max = max{A , ,4i}. This solu- 
tion, considered as a function of \ko) and |fci), will be 
denoted P^ ax (|fco),|fci)). 

The problem of unambiguously discriminating two 
non-orthogonal pure states \tpo), is a special case 
of the 1-dimensional kernel problem. Although cos 9 is 
defined to be the overlap of the 1-dimensional kernels, 
clearly KV'olV'i)! = cos9, and (JSJ becomes 



pmax = f 1 " 2^oPl1(^l)| if \{Mi>i)\ < yjf^ 
Pmax(l - \(ipo\i>i)\ 2 ) otherwise, 

_ (9) 
in agreement with [3|. In the case of equal prior proba- 
bilities, we have P max = 1 — \ (jpo\i)}i)\, as expected from 





IV. LOWER BOUND 

We consider a strategy that achieves UD of an arbi- 
trary pair of mixed states and which is strongly depen- 
dent on the geometrical relationship between the two sub- 
spaces /Co and /Ci : 

Theorem Consider two arbitrary mixed states po and 
pi . Denote the dimensionality of their kernels, /Co and 
JCi by ro and r\, and assume that ro > r\. There exist 
orthonormal bases for )Cb(b = 0,1) such that 

for 1 < j < r Q , l<i<n, 



{ki\k\) 



kj cos( 



(10) 
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where the Oj are the canonical angles between /Co ond /Ci 
}. In this case, the expression 

ri ro 

Pl=J2 P ™ X (\K)AH))+ E ( fc olPil^> (11) 



j=n+l 



forms a lower bound on the maximum probability of dis- 
criminating unambiguously between p Q and p\ . 

Proof The proof is constructive. Let Xf, be any nxrj 
dimensional matrix whose orthonormal columns span /Ct,. 
Define an rj, x rj, unitary matrix [/&, and an ro x n matrix 
S via a singular value decomposition 0] 



X\Xi = u sul 



(12) 



The matrix S is of the form S 

diagonal matrix of the form C = 
9i S [0,7r/2], while O is an (r - 



C 

o 

diag(cos9i, ... 
r\) x T\ matrix of O's 



where C is a 
cos 6 r 



Defining 



7t/2 for j > n, and denoting by |A$) the 



i'th column of X^Ub, we have constructed bases for /Co, 
JCi which satisfy (k J \kl) = SijeosOi, as required in the 
theorem. 

The measurement achieving the lower bound is associ- 
ated with a POVM wherein E a = TdU a i\ k \) i k \l E i = 
J2i=i Pi\K)(K\- The constraint Q takes the form 



fci)<fcll"E^l fc o><fco 



> 0. 



(13) 



and 



Since the 2 dimensional subspace spanned by 
\k{) for 1 < j < r\ is orthogonal to all other such sub- 
spaces and is orthogonal to \k J ) for r% < j < ro, the 
constraint H13|l reduces to constraints of the form (0) for 
1 < j < r i, and constraints of the form /— fjj \k 3 ) (fcp| > 0, 
for r\ < j < ro. In this manner, wc have reduced the 
problem to ri separate optimizations of the form already 
considered in section ITlI Bl Solving each of these yields 
the first term on the right hand side of I|11J) in the the- 
orem. The remaining ro — r\ optimizations are achieved 
by taking /3j = 1, which yields the second term on the 
right hand side of l|llfl . B 

In order to understand the geometry of the eigenbases 
for the POVM elements in this lower bound, it is helpful 
to realize that the canonical angles 9i form the unique ge- 
ometrical invariants describing the relationship between 
two subspaces. They can be defined iteratively: 6\ is the 
smallest angle between any pair of vectors drawn from 
/Co and K.±, and \kg), \k{) are the corresponding pair of 
vectors. 62 is the smallest such angle after these two vec- 
tors are removed, and so on. In this way, one obtains a 
simple geometrical picture of the measurement achieving 
the lower bound. We note that we have not found any 
example of UD wherein this lower bound is not optimal. 
Nonetheless, given that the eigenbases for E , Ex depend 
only on the subspaces spanned by po, p\ and not on the 
states themselves, there is no reason to expect it to be 
optimal in the general case. 



V. UPPER BOUND 

Theorem An upper bound on the maximum probabil- 
ity of unambiguously discriminating two mixed states, po 
and pi, is 



< 



{ 



l^^/poWi^Po, Pi) ifF(p ,pi) < 
Pmax(l - F(po, pi) 2 ) otherwise, 



(14) 

where F(po,pi) — Tr\ yfpo ^fp\\ is the fidelity. 

Proof Let be a purification of pt, [Tij . Clearly, UD 
of IV'o); l^i) with priors po,Pi can be achieved with a max- 
imal probability of success that is greater than or equal 
to the maximum probability of success in the UD of po,p\ 
with the same priors. This follows from the fact that any 
discrimination procedure for the latter task serves also 
as a discrimination procedure for the former task. Thus, 

P max ( Po po,PiPi) < , mm P max (p ^o)(^o|,Pi|^i)(^i|). 

where the minimization is over all purifications of pi,. 
We know, however, that P max (po|'0o)('0o| 5 Pi|V'i)(' i /'i|) i s 
given by Eq. (jSJ) . Minimizing either of the expressions on 
the right hand side of requires maximizing |('0o|'0i)l- 
Uhlmann's theorem 14] states that the maximum overlap 
between purifications of two density operators is equal to 
the fidelity between the density operators. Applying this 
observation yields the desired upper bound. ■ 

In the case of equal prior probabilities, the upper 
bound has the simple form p max < 1 — F(po,p±). Nu- 
merical studies of rank 2 mixed states in a 4 dimensional 
Hilbert space indicate that even for randomly chosen po 
and pi, our upper and lower bounds are generally very 
close. 



VI. SOME APPLICATIONS 
A. State comparison 

In 0, Barnett, Chefles and Jex introduced the follow- 
ing problem: given two systems, each of which is in one of 
two (generally non-orthogonal) states {\ipi} , \if) 2 }} , with 
what probability can one determine whether the two sys- 
tems are in the same state or in different states? Here, 
we consider the case when this determination must be 
made without error if at all. Assuming equal likelihood 
for either possibility, we can clearly interpret the problem 
as one of achieving UD of the two mixed states 

PO = ^IV'l'^lX^l^ll + ^|V , 2-02><-02^2| 

pi = lliPi^iiPiH + il^iPijifaH- (15) 

as was recognized by the authors of 0- 

These are rank two mixed states in a 4-dimensional 
space; as such their kernels are also two dimensional. 
We let \4>i) be the state orthogonal to \tpi), and choose 
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phases such that (^ll^) is real. It is clear that /Co 
is spanned by {IV'iV^j IV^i)} while /Ci is spanned by 
{\xp\_tpx) , \ip2ip2)}- Using the techniques of section (|IVJ| . 
we find that the canonical basis vectors in /Co are Ifcp) = 
I ^1^2) + \tp2ipi}_, \ko) =_feV>2) - |^ 2 -0i),_while_those of 
/Ci are = |^~i> + ^2^2), \kj) = |V>i^i>-|^2> (all 
suitably normalized). Since (fc 2 |fc 2 ) = 0, these two states 
can be included in their respective POVM elements with 
weight 1, i.e. 0:2 = 1 = 02- The remaining pair of states 
have overlap cos#i = (fcplfcj), and we can use the solu- 
tion for 1-dimcnsional kernels. Performing this calcula- 
tion, we find that our lower bound is 1 — F{p^ 1 p\). Since 
this is equal to our upper bound, we automatically know 
that this is optimal, and we have the optimal POVM by 
construction. 



B. State filtering 

A problem that has been considered recently by Sun, 
Bergou and Hillcry [pj is that of unambiguously dis- 
criminating whether a state is or whether it is in the 
set {|^2) , 1^3)} when the prior probabilities of the three 
states are 771 , 772, and rjs respectively. This task has been 
called state filtering. It is straightforward to see that it 
is an instance of UD of mixed states. Specifically, it is 
the problem of unambiguously discriminating 

po = \th) (H , pi = ^ (H + ^ m (H , 

with prior probabilities po — 771, pi = r\2 + V3- It is easy 
to find our lower bound for this problem. Since JCi is 
1-dimensional, we associate it with a unique vector \k\). 
The kernel /Co is 2-dimensional and so there is flexibil- 
ity in the basis that diagonalizes E\. Our lower bound 
dictates that we use the basis |&Q),|fcg), where \k^) is the 
vector in /Co that is maximally parallel to |fci), and \k^) is 
the vector in /Co that is orthogonal to both \k\) and \k^). 
Our lower bound is P L = P^{\kl),\k x )) + (fcg|pi|fcg) 
Defining the canonical angle 9 by cos 9 = (fci |/cq) and ap- 
pealing to the geometry of the problem, we find Aq = 
p (l — cos 2 9) and A± — piF(p , pi) 2 tan 2 9. Finally, 
defining F = F(po, pi)/ \/po/pi, the lower bound takes 
the following form: 

r l-p Q cos 2 9- ^F(p , Pl ) 2 F<cos 2 9 
P max > <^ 1 - 2^plF(p , pi)) cos 2 9 < F < 1 
I Pl {l-F{p , Pl ) 2 ) F>1 



This coincides precisely with the optimal solution derived 
in [lfj. Thus, our lower bound is found to be optimal in 
this case. 

By recasting state filtering as an instance of UD of a 
pure state and a mixed state, it has become apparent 
that the particular ensemble of states from which the 
mixed state is formed (in this case, ^2) an d with 
priors r\2 and 773) is not significant. As such our analysis 
applies also to the case wherein the second set of states 
contains an arbitrary number of elements, in agreement 
with the solution found in [llj . Moreover, the mixed 
state need not even arise as a result of ignorance of the 
pure state, but rather may arise as a result of the system 
being entangled with another degree of freedom. This 
is confirmed by the fact that the probability of success 
depends only on p\ . The problem of UD between multiple 
sets of pure states is a straightforward generalization 
of state filtering that can be recast as a problem of UD 
of mixed states in a similar way. 



VII. CONCLUSION 

We have shown that UD of a pair of mixed states is 
possible when these have distinct supports. We have pro- 
vided the maximum probability of success for the case 
of mixed states with orthogonal or 1-dimensional ker- 
nels. In the general case, we have determined an upper 
bound based on the probability of UD for purifications 
of the states, and a lower bound based on the geometri- 
cal invariants between their kernels. In addition to the 
tasks discussed here, there are likely to be many others 
for which UD of mixed states is useful, such as quantum 
random access codes, quantum oblivious transfer and en- 
tanglement distillation, to name a few. 
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